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ABSTRACT 


The  vertical  structure  of  the  mean  longshore  current  is  examined  on  three  strong 
current  days  during  the  DUCK94  experiment  and  it  is  found  well  described  by  a  logarithmic 
profile  (mean  correlation  coefficient  for  all  22  profiles,  0.98).  This  hypothesis  works  better 
in  the  trough  where  turbulent  bottom  boundary  layer  processes  are  predominant  than  over 
the  bar,  where  breaking-wave  induced  turbulence  generated  at  the  surface  modifies  the 
profile. 

A  relationship  between  near  bottom  vertical  velocity  profiles  and  bottom  roughness 
was  found.  The  bed  shear  stress  coefficient  varied  by  an  order  of  magnitude  across  the  surf 
zone  (0.001-0.05).  For  the  three  days  considered,  it  is  concluded  that  the  bed  shear  stress 
coefficient  increased  with  increasing  bottom  roughness,  and  therefore  is  an  important 
parameter  to  characterize  the  bottom  boundary  layer. 

The  influence  of  the  wind  is  parameterized  by  an  exponential  approximation  that 
works  well  in  50%  of  the  analyzed  profiles.  This  parameterization  assumes  that  the  residual 
Data-Normal  distribution  of  data  is  only  due  to  the  wind  and  does  not  consider  the 
alongshore  component  of  mass  transport  velocity,  which  is  assumed  small. 

These  results  should  be  treated  carefully  since  there  were  obtained  for  only  three  days 
over  a  barred  beach  with  strong  longshore  currents,  and  cannot  be  generalized  without 
further  studies. 


V 


TABLE  OF  CONTENTS 


l.  INTRODUCTION  . 1 

II.  THEORY . 5 

A.  BOTTOM  BOUNDARY  LAYER .  5 

B.  SURFACE  BOUNDARY  LAYER .  8 

m.  DUCK  94  EXPERIMENT  . 15 

IV.  DATA  RESULTS  . 19 

V.  DISCUSSION . 25 

VI.  CONCLUSIONS  . 29 

APPENDIX  A.  FIGURES . 31 

APPENDIX  B.  TABLES . 57 

LIST  OF  REFERENCES . 61 

INITIAL  DISTRIBUTION  LIST . 63 


VI 1 


LIST  OF  FIGURES 


1.  SLED  on  the  Beach  During  DUCK94  Experiment . 33 

2.  Bathymetry  and  Equipment  Positioning  During  DUCK94  Experiment  (1 1  Oct,). 

Sled  Transect  Line  was  940  m  Alongshore . 34 

3.  Bathymetry  Profiles  for  the  Three  Days  Considered  (10,  1 1  and  12  Oct.) . 35 

4.  Climatology  for  the  Three  Days  Considered  (10,  11  and  12  Oct  ).  Currents  Were 
Measured  in  the  Middle  of  the  Trough.  H  is  the  Significant  Wave  Height  and  T  is 

the  Period  of  Peak  Frequency . 36 

5.  (a)  Least-Squares  Fit  to  the  Data  for  10  Oct . 37 

5.  (b)  Least-Squares  Fit  to  the  Data  for  11  Oct . 38 

5.  (c)  Least-Squares  Fit  to  the  Data  for  12  Oct .  39 

6.  Example  of  Vertical  Profile  (the  Upper  Two  Current  Meters  Sometimes  Were  Out 

of  the  Water) . . 40 

7.  (a)  Vertical  Profiles  of  Mean  Longshore  Currents  for  10  Oct . 41 

7.  (b)  Vertical  Profiles  of  Mean  Longshore  Currents  for  1 1  Oct . 42 

7.  (c)  Vertical  Profiles  of  Mean  Longshore  Currents  for  12  Oct . 43 

8.  Variation  of  the  Rms  Bottom  Roughness  with  the  Cross-Shore  Distance . 44 

9.  Relationship  Between  c ^  and  /  V . 45 

z 

10.  Relationship  Between  and  /  V . 46 


IX 


11.  Bed  Shear  Stress  versus  Bottom  Roughness.  The  Line  Represents  a  Linear 

Regression  Based  on  All  Points  but  the  Outlier .  47 

12.  c,  versus  — . 48 

^  A 

z  j- 

13.  — ^versus  —  (Outlier  not  Included) . 49 

h  A 

14.  versus  Cross-Shore  Distance . 50 

15.  V  ^  versus  Cross-Shore  Distance .  51 

16.  Relationship  Between  andy . 52 

17.  Relationship  Between  v^.  andy . 53 

18.  Example  of  Wind  Contribution  to  the  Longshore  Current:  Exponential 

Approximation  to  the  Residual .  54 

z  Hr 

19.  Relationship  Between — -  and  — ^  — . 55 

h  h 


X 


LIST  OF  TABLES 


1.  Profile  Fitting  Results . 59 

2.  Wind  Data . 60 


XI 


ACKNOWLEDGMENTS 


The  author  wishes  to  express  his  appreciation  to  all  those  who  participated  in 
DUCK94  experiment,  namely  the  staff  of  the  U.S.  Army  Field  Research  Facility  under 
the  direction  of  B.  Birkemeier.  In  addition,  special  appreciation  is  expressed  to  R.  Wyland, 
Naval  Postgraduate  School  for  his  role  in  acquisition  of  wave  and  current  data,  and  to 
Mary  Bristow,  Naval  Postgraduate  School,  for  help  in  processing  data.  Thanks  to  LT  Jeff 
Swayne,  USN,  for  providing  bottom  roughness  data.  Finally  a  special  thanks  to  Tom 
Lippmann,  Naval  Postgraduate  School,  to  Timothy  Stanton  and  Edward  Thornton  for  their 
important  suggestions  and  orientations. 

Lastly  and  most  importantly,  I  want  to  express  my  sincere  thanks  to  my  wife,  Sonia. 
Without  her  love,  support,  encouragement  and  understanding,  this  thesis  would  never  have 
been  written. 


xiii 


I.  INTRODUCTION 


Knowledge  of  the  bottom  and  surface  boundary  layers  is  fundamental  to 
understanding  nearshore  hydrodynamic  and  sediment  processes.  For  steady  flow,  such  as 
in  a  river,  the  bottom  boundary  layer  is  well  described  by  a  logarithmic  profile.  The 
superposition  of  waves  on  the  mean  current  produces  enhanced  bottom  friction  (e.g..  Grant 
and  Madsen  (1979),  Christofferson  and  Jonsson  (1985),  Myrhaug  and  Slaattelid  (1989)  and 
Sleath  (1990  ) ).  As  a  result,  the  gradient  of  mean  current  near  the  bed  is  decreased  and  a 
more  linear  profile  can  be  expected.  The  influence  by  waves  is  inversely  proportional  to  the 
depth  of  the  water  with  decreasing  importance  of  wave-bottom  boundary  with  increasing 
depth. 

Only  limited  wave-driven,  longshore  current  vertical  profiles  have  been  measured. 
Visser  (1984)  measured  wave-driven  longshore  currents  in  a  laboratory  experiment  using 
dye  displacement  over  the  vertical.  He  found  the  alongshore  spatially  averaged  currents 
showed  more  linear  than  logarithmic  profiles.  However,  Visser  (1986),  repeated  the 
experiment  in  the  same  basin,  measured  the  vertical  profiles  using  micro-propeller  and  laser 
Doppler  anemometers  and  found  profiles  approached  more  a  logarithmic  form.  In  a  similar 
laboratory  wave-driven  longshore  current  experiment  using  laser  anemometers,  Simons  et 
a/.  (1992)  found  the  vertical  velocity  profiles  tend  to  show  logarithmic  behavior.  Church  et 
al.(\993)  used  field  data  from  SUPERDUCK  experiment  (1986)  to  study  the  problem.  Three 
electromagnetic  current  meters  mounted  on  a  mobile  sled  were  used  in  a  vertical  stack  0.7, 
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1.0  and  1.5  m  above  the  bed.  The  results  were  not  conclusive  in  relation  to  the  logarithmic 
profile  question  since  only  three  current  meters  were  used.  It  was  concluded  that  more 
current  meters  located  nearer  the  bed  are  required  to  determine  the  bottom  boundary  layer 
shape. 

Turbulence  induced  by  breaking  waves  can  modify  the  vertical  profile  of  longshore 
currents.  The  downward  vertical  momentum  mixing  produced  from  the  breaking  wave- 
injected  turbulence  results  in  a  more  uniform  velocity  profile,  increasing  the  shear  stress  at 
the  bed  (Deigaard  et  al,  1986  and  Church  and  Thornton,  1993). 

Wind  also  modifies  the  vertical  profile  of  longshore  currents.  Murray  (1975) 
derived  the  vertical  profiles  of  mean  currents  as  a  function  of  the  wind  in  shallow  water  by 
balancing  the  surface  stress  by  wind  with  the  bottom  stress,  with  and  without  the  Coriolis 
force.  He  compared  the  theory  with  observations  of  current  profiles  obtained  using  drogues 
at  various  elevations  well  outside  from  surf  zone.  Comparing  wind  and  wave  forcing  of 
longshore  currents  within  the  surf  zone,  Whitford  and  Thornton  (1993),  showed  the  wind 
can  be  significant,  but  only  at  high  speed  and  when  the  wind  has  a  strong  alongshore 
component. 

In  the  following,  mean  longshore  current  profiles  obtained  over  a  barred  beach  are 
examined  with  the  following  objectives: 

a.  Test  the  hypothesis  that  the  turbulent  bottom  boundary  layer  of  the  mean  longshore 
currents  is  logarithmic. 

b.  Determine  how  breaking- wave  induced  turbulence  generated  at  the  surface 
modifies  the  vertical  velocity  profile. 
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c.  Investigate  how  wind  modifies  the  vertical  profile  of  the  longshore  currents. 

d.  Examine  the  relationship  between  near  bottom  vertical  velocity  profiles  and 
bottom  roughness,  including  the  influence  of  ripples  and  mega-ripples. 
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II.  THEORY 


The  vertical  profile  of  the  longshore  currents  is  significantly  affected  by  the  physical 
boundary  layers  of  the  bottom  and  sea  surface.  The  bottom  boundary  assumes  more 
importance  since  it  determines  the  general  logarithmic  profile  shape.  However,  processes 
in  the  surface  layer  can  modify  the  profile  in  the  presence  of  waves  and  wind. 


A.  BOTTOM  BOUNDARY  LAYER 

Neglecting  molecular  viscous  stresses,  the  alongshore  momentum  equation  (  y- 
direction)  is  written 


dpv  dpuv  8pv^  8pwv  dp  ,  . 

dt  dx  dy  dz  dy 


Let  M.=  U.+w  .+«.  and  m'=  w'  +w.,  where  the  horizontal  («.,  /=1,2)  and  vertical 

lilt  I  l’  \  |5  5  / 


w 


velocities  have  been  expanded  into  mean,  turbulent  and  wave-induced  components.  Assume 

0  _ 

straight  and  parallel  contours,  i.e.,  — (  )=0  (overbar  indicates  time  averaging),  neglect 

horizontal  turbulent  and  mean  momentum  mixing  ^  and  neglect  since  its 

ax  8x  dz 

contribution  is  small  above  the  wave  bottom  boundary  layer.  Assuming  steady  state 


conditions  and  after  time  averaging 


dpvu  _  _dpw'v' 
dx  dz 


(2) 
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which  says  the  cross-shore  change  in  the  wave  induced  alongshore  momentum  flux  is 
balanced  by  vertical  changes  in  alongshore  turbulent  shear  stress.  In  the  shallow  water  of  the 
surf  zone,vz7  can  be  assumed  independent  of  depth,  and  its  gradient  in  x  is  constant  in  the 
>’-direction  for  straight  and  parallel  contours.  Defining 


t 


-piv  V  =  pVj 


dz 


(3) 


where  is  the  turbulent  eddy  viscosity.  Substituting  (3)  into  (2) 

=  2? 

dz  dx 


(4) 


where  is  constant  alongshore.  Constant  shear  stress  distributions  occur  in  flows  driven  by 
constant  hydrostatic  pressure  gradients,  such  as  in  steady  open  channel  flow,  which  are  well 
described  by  a  logarithmic  velocity  profile.  Therefore  it  is  hypothesized  that  the  vertical 
velocity  profile  for  longshore  currents  is  logarithmic.  A  steady,  uniform,  turbulent  boundary 
layer  in  the  alongshore  direction  can  be  described  by  the  logarithmic  “law  of  the  wall”: 

v(z)  =  —  In  ^  (5) 

K  Z 

£> 

where  z  is  positively  upward  from  the  surface,  h  is  the  depth  of  water,  k  is  the  Von 
Karmann  constant  (0.4),  is  the  shear  stress  velocity  and  z^  is  the  roughness  height.  The 
effect  of  including  waves  on  the  bottom  stress  is  to  effectively  increase  the  apparent 
roughness  height  denoted  by  replacing  z^  with  (that  is  analogous  to,  but  larger  thanz^ 
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since  it  includes  the  additional  stress  caused  by  the  waves).  Following  the  empirical  model 
by  Nielsen  (1992)  for  current  profiles  in  the  presence  of  waves,  can  be  related  to  the 
thickness  of  wave-dominated  layer  6  by 

z^  =  e  '^6  (6) 


Since  6  has  no  generalized  formulation,  a  particular  formula  is  suggested.  For  fairly  rough 

conditions  (0.06  <  —  <0.5  ),  constant  eddy  viscosity,  and  neglecting  the  angle  between 
A 

wave  and  current’  8  can  be  parameterized  as 


6  » 


0.016  Ao3  ^ 


K 


♦ 

V 


0) 


where  o)  is  the  radian  frequency,  r  is  the  hydraulic  roughness  and  A  is  the  orbital  amplitude 
of  the  waves  just  above  the  bottom  wave  boundary  layer. 

The  shear  stress  profile  is  determined  by  integrating  (4)  and  assuming  the  stress  at 
the  surface  is  equal  to  zero,  to  give 

f  (8) 

h 

where  h  is  the  total  depth.  Bottom  stress  is  related  to  v,  through 

(9) 


Substituting  (5)  and  (8)  into  (3)  we  obtain 
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(10) 


V,  =  -tv,r(l+i) 
h 

which  says  the  ediJy  viscosity  is  parabolic  with  depth.  However,  since  it  is  based  on  the  “law 
of  the  wall”,  (10)  would  not  be  expected  to  be  an  adequate  formulation  near  the  surface, 
where  the  wave  and  wind  effects  can  be  significant. 

B.  SURFACE  BOUNDARY  LAYER 

The  surface  boundary  layer  is  modified  by  the  waves  and  wind.  Wave  modifications 
can  be  due  to  three  contributions:  modification  of  mean  currents  by  measuring  currents  in 
a  Eulerian  frame  with  an  undulating  boundary;  contribution  to  longshore  currents  of  mass 
transport  velocity  in  the  crest-trough  region  by  obliquely  incident  waves  in  a  Eulerian  frame, 
and  modification  by  the  wave  breaking. 

The  most  significant  modification  of  the  near  surface  mean  current  profile  is  due  to 
the  effect  of  measuring  currents  in  an  undulatory  wave  boundary  in  a  Eulerian  frame.  In  the 
wave  crest/trough  region,  the  averaged  measured  velocity  is  decreased  because  the  current 
meter  is  out  of  the  water  part  of  the  time.  For  instance,  the  current  meter  is  out  of  the  water 
half  the  time  at  the  MSL,  and  the  time  averaged  current  is  only  50%  the  expected  value  of 
the  logarithmic  profile.  To  account  for  this,  the  surface  elevation  probability  distribution 
function  (pdf)  is  applied  to  the  expected  mean  current  profiles  in  the  absence  of  waves.  The 
percent  of  time  the  current  meter  is  in  the  water  is  given  by  1-P(t)),  where  P(r|)  is  the 
cumulative  surface  elevation  pdf.  In  a  Eulerian  frame  of  reference,  the  modified  mean 
current  in  the  crest-trough  region  is  given  by 
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V(z)  =  [l-P(il)]V(z) 


(11) 


where  V(z)  represents  the  log  profile  in  the  absence  of  waves.  For  mild  wave  heights  and 
in  deeper  water,  the  surface  elevation  pdf  is  well  described  by  the  Gaussian  pdf.  As  will  be 
seen,  the  measured  pdfs  in  the  surf  zone  are  slightly  positively  skewed  from  the  Gaussian 
distribution. 

The  mass  transport  velocity  can  contribute  to  the  longshore  current  in  the  upper 
boundary  layer  when  the  incident  wave  angle  is  relatively  large.  Assuming  irrotational  flow 
and  evaluating  mass  flux  (transport),  M,  by  considering  separately  two  regions  in  a  Eulerian 
frame  of  reference  (cf.  Philips,  1977) 

■^  =  P  +  Qudz)  (12) 

The  first  term  is  zero  for  irrotational  flow.  For  infinitesimal  linear  wave  theory,  u  is  not 
defined  between  zero  and  the  surface  r| .  A  Taylor  series  expansion  about  z=0  is  used  to 
extend  defined  values  of  u(0)  to  the  surface,  giving 

M  =  p«(0)Ti  (13) 

correct  to  0(a/k)^ .  Applying  linear  theory 

M  =  A  04) 

pC 

where  E  is  the  energy  and  C  the  celerity.  The  mass  transport  for  a  single  wave  is  interpreted, 
to  this  order  in  an  Eulerian  reference  frame,  as  due  to  a  uniform  velocity  confined  to  the 
crest/trough  region.  The  mass  transport  velocity  is  defined 
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u 


M 

plA 


(15) 


where  A  is  the  amplitude. 

For  random  waves,  the  wave  amplitudes  are  reasonably  described  by  the  Rayleigh 
distribution.  The  only  waves  that  contribute  at  any  elevation  2,  will  have  an  amplitude 
A  >  z.  Hence,  the  ensemble  average  mass  transport  velocity  profile  is  obtained  by  applying 
the  wave  amplitude  pdf 


<U(z)>  =  r  U(A)p(A)dA 
Jz^A 


(16) 


The  Rayleigh  distribution  is  given  by 


P(.A)  = 


ZA 

H  2 


JA  2 

rmj 


(17) 


where  is  the  rms  wave  height.  Substituting  (17),  U(z)  is  given  by 


<Z7(z)> 


CO  r  ^ 

- [ -  e 

8  tanh  kh  H 

rms 


H- 


2z 


rms 


■)] 


(18) 


where  erfc(x)  is  the  complementary  error  function.  The  alongshore  component  is  defined  by 

<V(z)>  =  <U(z)>sinQ  (19) 


where  0  is  the  mean  wave  incident  angle. 
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The  surface  boundary  layer  during  wave  breaking  is  typified  by  an  intense 
production  of  turbulence  that  eventually  is  dissipated  in  the  shear  layer  at  the  lower 
boundary  of  the  surface  roller.  The  remaining  energy  propagates  vertically  through  the  water 
column.  The  increase  of  eddy  viscosity  due  to  breaking-induced  turbulence  produces  a  more 
vertical  profile  of  the  current  for  a  given  shear  stress  distribution,  compared  with  profiles 
in  the  absence  of  breaking.  This  results  in  a  stronger  velocity  gradient  near  the  bed, 
increasing  the  bed  shear  stress. 

The  influence  by  wind  on  the  longshore  currents  can  be  significant  during  times  of 
strong  alongshore  wind.  Again  integrating  (4),  but  including  an  alongshore  wind  stress 

^  (11 -^)  (20) 


withi:^=  (surface)  and  T^=  pC^|ir|vj_  (bottom),  where  W  is  the  wind  speed, 

is  the  density  of  air,  (f)  the  wind  angle  with  the  longshore  current,  the  drag  coefficient, 
the  bed  shear  stress  coefficient  andvj_  the  longshore  current  speed  at  the  bottom.  Equating 
the  stress  at  the  surface  with  the  alongshore  wind  stress  component  we  have 


(21) 


Solving  for  v  exclusively  due  to  the  wind 


V 


-^C^Fr"cos(t)j“ 

P 


(22) 
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Murray  (1975)  considered  constant 


P 


(23) 


This  theory  applies  when  the  wind  angle  to  the  coast  is  less  than  45  degrees  and  does  not 
consider  the  influence  of  the  waves.  It  is  found  this  treatment  overestimates  the  speed  at  the 
bottom,  which  was  verified  using  field  data  from  DUCK94. 

More  complex  forms  can  be  adopted  for  eddy  viscosity  such  as  the  half-parabolic 
distribution  (Roelvink  and  Reniers,  1994).  Herev^  is  defined  as 


1  ^  i  (2-i) 

2  ‘  h  h 


with 


As  will  be  seen,  an  exponential  distribution  of  the  wind-driven  velocity  appears  to 
give  a  reasonable,  fit  to  the  data 

V  =  (26) 

where  is  the  velocity  at  MSL  (top  current  meter)  and  P  a  free  parameter.  The 
corresponding  eddy  viscosity  is 
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(27) 


^  P(z-*) 

P  V,  P 


This  approach  is  used  to  interpret  the  data  from  DUCK94. 
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III.  DUCK  94  EXPERIMENT 


The  measurements  described  here  are  part  of  the  comprehensive  nearshore  DUCK94 

experiment  conducted  during  August  and  October  1994  at  the  U.S.  Army  Corps  of 
Engineers  Field  Research  Facility  (FRF),  Duck,  North  Carolina.  The  FRF  is  located  on  the 
Outer  Banks,  a  barrier  island  formation  with  no  major  coastal  structures  that  can  obstruct 
nearshore  flows.  The  beach  is  a  two-bar  system  with  a  dynamic  inner  bar  (30-120m 
offshore)  and  a  secondary  bar  with  lower  amplitude  (300-400m  offshore).  The  mean 
foreshore  slope  of  the  beach  is  approximately  0.08  (1 : 12)  and  the  slope  offshore  of  the  bars 
is  approximately  0.006  (1:170)  (Lippmann,  1993).  The  mean  tidal  range  is  1.0  m.  Sediments 
within  the  surface  zone  are  well  sorted  with  a  mean  grain  size  of  0.2  mm.  Sediments  on  the 
foreshore  are  poorly  sorted  with  larger  mean  grain  size  (>0.4  mm). 

A  specially  designed  sled  was  used  as  a  platform  to  mount  instruments  (Fig.  1).  The 
sled  is  constructed  of  a  3  x  4  m  six-inch  aluminum-pipe  frame  with  two  5m,  8  inch  pipe 
runners.  This  low-profile  structure  was  stabilized  by  180  Kg  of  lead  weight  plus 
approximately  450  Kg  of  sand  inside  the  runners.  In  addition,  four  fins  (45  cm  wide) 
extended  60  cm  into  the  sand  to  insure  the  sled  did  not  move  in  the  surf  zone. 

The  vertical  profile  of  longshore  current  was  measured  using  a  vertical  stack  of  8 
Marsh-McBimey  two  component  electromagnetic  current  meters  with  2.5  cm  diameter 
spherical  probes  mounted  on  a  2.5  m  mast  (Fig.  1).  The  current  meter  elevations  above  the 
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bed  were  23,  42,  68,  101,  147,  179,  224  and  257  cm.  The  sled  was  oriented  such  that  the 


vertical  slack  of  current  meters  was  consistently  placed  on  the  “up-current”  side  of  the  sled 
to  avoid  flow  contamination  by  the  sled  structure  during  observations.  The  current  meters 
were  pre-  and  post-  calibrated  in  a  tow  tank  at  the  Naval  Postgraduate  School  with  an 
agreement  of  1.9  %  in  gain.  An  in  situ  determination  for  the  offset  is  used,  which  was 
obtained  by  reversing  the  orientation  of  the  current  meters  on  a  very  slow  longshore  current 
day  (8  October)  by  turning  the  sled  around  and  returning  it  to  the  same  location  (within  1 
m)  within  one  hour.  The  in  situ  determined  offsets  were  within  1  cm/s. 

Waves  and  mean  water  level  were  measured  using  an  array  of  six  pressure  sensors 
configured  in  a  3m  square  with  one  in  the  center  and  another  in  the  middle  of  a  cross-shore 
leg. 


The  data  Were  digitally  encoded  to  14  bit  precision  on  the  sled  at  36  samples/second 
and  transmitted  to  shore  via  a  fibre-optic  cable  where  signals  where  monitored  and  recorded. 
Short  cables  from  the  sensors  to  the  data  acquisition  system  on  the  sled  (<7m)  resulted  in 
exceptionally  low  noise  electromagnetic  current  meter  and  pressure  sensor  signals.  An 
armored  cable,  married  to  the  sled  chain  tether,  provided  power  and  controller  signals  for 
the  instruments  via  two  conductors  and  returned  the  digitized  signals  and  video  via  two 
fibre-optic  lines. 

The  sled  was  towed  to  the  farthest  offshore  location  for  the  first  run  (approximately 
160m  from  the  shoreline)  by  the  Coastal  Research  Amphibious  Buggy  (CRAB).  A  four- 
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wheel  drive  forklift  pulled  the  sled  shoreward  10  to  30  meters  for  subsequent  runs  (each  run 
was  nominally  one  hour).  Five  to  eight  runs  were  made  across  a  transect  each  day.  The  data 
were  acquired  during  daylight  to  early  night. 

The  morphology  of  the  bottom  (bathymetry)  was  measured  at  various  scales  from 
the  CRAB.  Large-scale  variations  of  bathymetry  were  obtained  by  using  an  autotracking 
laser  ranging  system  to  measure  the  CRAB  position  with  a  vertical  accuracy  of  less  than  3 
cm.  Contour  plots  of  the  bathymetry  for  the  days  selected  for  analysis  show  the  contours  in 
the  vicinity  of  the  sled  measurements  to  be  quasi  straight  and  parallel  (Fig.  2;  October  1 1th 
is  typical).  SmalLscale  vertical  bottom  variations  relative  to  the  CRAB,  including  ripples 
and  megaripples,  were  acquired  with  a  1  MHz  sonic  altimeter  mounted  on  the  CRAB  70  cm 
from  the  bed  with  a  vertical  accuracy  less  than  2  cm  (Gallengher  et  al,  1995).  Bathymetry 
for  selected  days  (Fig.  3)  shows  a  pronounced  bar  progressively  moving  offshore  and 
significant  small-scale  morphology  in  the  trough.  Areal  variations  were  determined  using 
a  500  Khz  side-scan  sonar  also  mounted  on  the  CRAB. 

Meteorological  information  of  wind,  air  temperature,  atmospheric  pressure  and  sea 
surface  temperature  were  recorded  simultaneously  at  the  end  of  the  FRF  pier  and  onshore 
near  the  FRF  laboratory. 
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IV.  DATA  RESULTS 


The  data  analysis  presented  is  for  the  October  phase  of  the  DUCK94  experiment.  The 
weather  during  October  is  climatologically  characterized  by  three  distinct  phases  (Fig.4): 
weak  currents  and  winds  from  north  (4-9  Oct),  relatively  strong  currents  from  north  (0.6- 1.0 
m/s)  caused  by  a  storm  with  predominant  winds  and  waves  from  north  (10-12  Oct.),  and 
variable  currents  and  winds  from  north/south  (13-21  Oct.).  The  data  selected  to  analyze  are 
10-12  October,  during  the  strong  longshore  currents  period. 

All  22  vertical  profiles  of  longshore  currents  obtained  during  the  three  days  are  used. 
The  profiles  are  based  on  the  measurements  by  seven  current  meters  over  the  vertical.  The 
current  meter  the  closest  to  the  sea  bed  was  not  used  because  of  malfunction.  Mean 
longshore  velocities  are  1-hour  averaged  data  with  the  exception  of  three  runs,  run  7  on  1 1 
Oct.  and  runs  6  and  7  on  12  Oct.,  which  are  40-minutes  averaged  data.  In  examining  the 
bottom  and  surface  boundary  layers,  the  current  meter  data  are  treated  separately.  For  the 
bottom  boundary  layer,  logarithmic  profiles  based  on  a  least  square  method  were  fit  to  the 
data  to  test  the  validity  of  the  logarithmic  profile  hypothesis  (Fig.5).  The  linear  correlation 
coefficients  ranged  from  0.89  to  0.99  (Table  1).  Since  the  logarithmic  velocity  profile 
hypothesis  is  a  bottom  boundary  layer  concept,  information  from  current  meters  near  the 
surface  can  be  contaminated  by  the  influences  of  waves  and  wind  and  coming  in  and  out  of 
the  water  and,  hence  were  not  included.  A  criterion  is  used  such  that  only  information  from 
current  meters  below  (MSL-  H^^),  where  MSL  stands  for  mean  sea  level,  is  considered  to 
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define  logarithmic  profiles  of  the  bottom  boundary  layer.  This  criterion  assures  the  current 
meters  in  the  analysis  came  out  of  the  water  less  than  0.25%  of  the  time  based  on  a  Gaussian 
distribution,  which  is  conservative  for  the  measured  skewed  distributions.  The  rms  wave 
height  is  approximated  by  where  is  the  variance  calculated  from  the  surface 

elevation  time  series.  The  surface  elevation  was  calculated  by  Fourier  transforming  the  one- 
hour  pressure  record,  applying  a  linear  wave  theory  transfer  function  to  the  complex  Fourier 
amplitudes  in  the  frequency  domain,  and  inverse  transforming  to  obtain  the  surface  elevation 
time  series. 

In  the  surface  boundary  layer,  the  upper  current  meters  sometimes  came  in  and  out 
of  water,  which  causes  noise  in  the  current  meter  output.  To  eliminate  this  noise,  when  the 
current  meters  were  within  5  cm  beneath  the  surface  or  higher  as  determined  from  the 
pressure  sensor  time  series,  the  current  velocities  were  set  to  zero. 

An  example  of  combining  the  surface  and  bottom  boundary  layers  is  given  in  Fig. 
6.  The  bottom  is  described  by  a  logarithmic  profile  and  the  surface  boundary  layer  by  the 
measured  surface  elevation  pdf  (equation  11). 

The  velocity  profiles  at  successive  offshore  positions  (runs)  that  the  sled  occupied 
during  a  transect  are  shown  in  Fig.  7.  The  data  agree  well  with  the  logarithmic  profile.  The 
largest  discrepancies  occur  over  the  bar,  where  wave  breaking  is  most  intense  (x-distance 
between  220  and  240  m).  The  largest  deviation  occurs  for  the  fourth  run  of  1 1  Oct,  which 
generated  an  outlier  in  the  parameter  estimates.  The  average  correlation  for  all  u  profiles 
between  log  elevation  and  longshore  velocity  is  0.98;  this  value  is  decreased  by  the 
observations  in  the  breaking  zone  (over  the  bar). 
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The  value  of  is  calculated  from  the  intersection  of  the  linear  regression  least 
squares  fit  to  the  data  in  semi-log  form,  and  the  shear  stress  velocity  is  estimated  from 
the  slope.  In  addition,  the  bed  shear  stress  coefficient,  C^,  can  be  calculated  from  the  data. 
Assuming  a  quadratic  bed  shear  stress  relationship 


=  pC^(«2+  v^)^v 


(28) 


and  combining  with  (8),  gives 

2 

c  =  _ _  (29) 

"  (u^.vY^v 

A  mean  eddy  viscosity  for  each  run  is  defined  by  integrating  (10) 

V,  =  A  (30) 

6 

Values  of  v^,  andv^  are  given  in  table  1. 

An  attempt  is  made  to  find  relationships  between  the  calculated  parameters  C^,  z^ 

and  with  physical  measurements  such  as  (ratio  of  near-bottom  wave  velocity 

magnitude  to  mean  current  measured  at  lowest  current  meter),  y  =  (%  of  wave 

h 

breaking)  and  the  rms  bottom  roughness,  r,  measured  by  the  CRAB’s  altimeter.  An  example 
of  the  variation  of  rms  bottom  roughness  calculated  over  20  m  detrended  sections,  in  one 
meter  steps  across  the  surf  zone  is  shown  in  Fig.  8.  Bottom  roughness  is  smoothest  offshore 
and  over  the  bar  with  increased  roughness  within  the  trough. 
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The  data  are  qualitatively  sorted  by  location  over  the  bar,  in  the  trough  and  on  the 
foreshore.  This  method  allows  a  better  identification  of  the  possible  correlations  among 
variables. 

In  general,  and  values  increased  an  order  of  magnitude  across  the  surf  zone. 

U 

values  were  lowest  over  the  bar  and  tended  to  vary  over  a  small  range  for  all  values  of  — 

V 

.  .  .  w  ^ 

(Fig.  9).  tended  to  increase  with  increasing  —  over  the  trough  area.  —  tended  to 
^  V  h 

behave  in  a  similar  manner  to  (Fig.  10).  does  increase  with  bottom  roughness 

V 

(Fig.  11).  C.  is  also  plotted  against —  (Fig.  12)  where  A  is  the  wave  orbital  excursion 

A 

amplitude,  which,  for  linear  shallow  water  wave  theory  is  given  by 


where  g  is  the  gravity  constant.  —  also  tends  to  increase  with  —(Fig.  13). 

h  A 

Analyzing  the  parameters  over  the  cross-shore  distance,  increases  slightly  toward 

the  beach  but  this  smooth  trend  is  interrupted  with  a  sharp  peak  in  the  wave  breaking  area. 

After  this  returns  to  the  previous  behavior  (Fig.  14).  shows  minimum  values  over  the 

bar,  a  peak  in  the  trough  region  and  decreasing  values  on  the  foreshore,  reflecting  that  it  is 

Z  H  fTtt  s 

calculated  as  the  product  of  v^/i(Fig.l5).  Plots  of  C,  or  —  versus  y-  - show  poor 

^  h  h 

correlation  (Figs.  16  and  17). 

The  wind  contribution  to  the  longshore  currents  was  evaluated  by  plotting  an 
exponential  approximation  (26),  used  to  model  the  wind-driven  portion  of  the  profiles,  and 
compared  with  the  residual  of  the  measured  surface  elevation  pdf  (Fig.  1 8).  The  exponential 
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approximation  shows,  in  general,  good  agreement  with  the  residual.  The  wind  data  is 
summarized  in  Table  2. 
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V.  DISCUSSION 


There  has  been  a  lack  of  field  data  to  describe  the  mean  vertical  velocity  profile  of 
longshore  currents.  Three  days  of  DUCK94  experiment  when  strong  longshore  currents 
occured  were  examined.  Mean  currents  obtained  using  seven  current  meters  spaced  from  41 
to  257  cm  above  the  sea  bottom,  were  well  described  by  a  logarithmic  profile.  The  validity 
of  the  logarithmic  approach  is  reflected  in  a  high  correlation  mean  value  of  0.98  for  all  22 
profiles,  between  the  log  elevation  of  the  current  meters  and  the  measured  velocities. The 
largest  deviations  occured  over  the  bar  and  the  foreshore,  related  to  the  increased  turbulence 
caused  by  wave  breaking. 

A  sensitivity  test  was  performed  to  determine  if  the  correlation  coefficient  changed 
with  the  number  of  current  meters  used  to  define  the  profile.  Using  either  all  the  current 
meters  available  under  the  criteria  established  in  chapter  4,  or  only  the  lowest  three  or  four, 
the  correlations  are  similar  with  maximum  changes  of  only  0.02.  This  demonstrates  the 
validity  of  the  logarithmic  approach  over  most  of  the  water  colunm.  In  the  same  way,  the 
values  of  do  not  significantly  change  when  calculated  with  only  the  lowest  current 
meters. 

In  the  breaking  zone  the  values  of  the  bed  shear  stress  coefficient  (or  bottom  friction 
coefficient),  are  in  the  order  of  1 0  while  the  rest  of  the  transect  shows  values  in  the 
order  of  10'^  (consistent  with  previous  experiences  in  the  same  beach). 
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There  is  some  uncertainty  of  the  exact  distance  of  the  current  meters  from  the  bed, 
because  the  sled  runners  sunk  into  the  sand  an  unknown  amount  depending  on  bearing 
capacity  of  the  bed  (theory  and  divers  suggest  3-8  cm)  and  the  undulating  bottom, 
particularly  over  mega-ripples.  To  test  the  sensitivity  of  due  to  the  uncertainty  in  depth, 
the  elevation  of  the  current  meters  was  shifted  ±20  cm  in  steps  of  1  cm  and  recalculated. 
An  average  value  of  aC^  =  O.OOl/lOcm  was  calculated  (5  to  30  %  of  calculated  C^.  The 
changes  of  are  less  significant  for  increased  depth  (-20  cm)  than  for  decreased  depth 
(±20cm). 

z  u 

and  —  show  surprisingly  poor  correlation  with  the  wave/current  parameter  — , 

z  u 

since  it  was  expected  an  increasing  of  C.  and—  with  increasing  values  of  —  (increasing 

^  h  V 

of  the  wave  influence  over  the  vertical  profile).  Equation  (6)  suggests  a  relationship  between 
z^  and  wave  parameters  .  For  shallow  water  waves  (6)  can  be  rewritten 


-1O.OI6 


/— 

/  rms 


V. 


(32) 


z 

Equation  (32)  combines  wave,  current  and  roughness  variables.  —  is  plotted  against 

h 

equation  (32)  expression  in  Fig.  19.  However,  including  wave  and  current  variables  with 
roughness  makes  the  correlation  withz^  worse. 

The  strongest  correlation  was  found  between  and  the  bottom  roughness,  with 
increasing  with  the  bottom  roughness.  Theoretically  this  is  expected  as  larger  roughness 
implies  larger  bottom  stress  due  to  form  drag,  and  consequently  larger  values,  z^  also 
increased  with  the  bottom  roughness  ( — ),  but  no  as  strongly  as  C..  The  C.  values  varied 

A  J  J 
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by  an  order  of  magnitude  across  the  surf  zone.  This  implies  that  the  bottom  roughness  has 
an  important  role  in  the  characterization  of  the  bottom  boundary  layer  by  affecting  the  value 
of  C^.  Roughness  values  increase  toward  the  beach  and  are  larger  in  the  trough.  The 
increased  roughness  of  the  bottom  in  the  trough  is  associated  with  mega-ripples  observed 
by  side-scan  sonar  (Swayne,  1995). 

The  behavior  of  v^  along  the  cross-shore  line  is  basically  a  function  of  the  velocities 
and  the  depth.  Low  values  over  the  bar  and  foreshore  (smaller  h  )  and  high  values  over  the 
trough  (larger  velocities)  are  explained  by  (30). 

Wind  contribution  to  the  longshore  current  was  examined.  Using  a  simplified 
approach,  reasonable  agreement  was  found  between  an  exponential  approximation  and  the 
residual  of  velocity/surface  elevation  pdf  correction.  Nevertheless  it  is  difficult  to  predict  if 
the  residual  is  due  only  to  the  wind  or  if  other  factors,  such  as  mass  transport  velocity  or 
wave  breaking  also  contribute  to  this  residual.  The  free  parameter  P  in  equation  (26)  varied 
between  3  and  4  giving  reasonable  predictions  for  the  cases  studied.  Since  the  approximation 
was  tested  over  only  12  profiles,  future  work  is  needed  to  refine  the  value  of  P . 
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VI.  CONCLUSIONS 


The  vertical  structure  of  mean  longshore  currents  on  a  barred  beach  is  well  described 
by  a  logarithmic  profile  for  the  three  strong  longshore  current  days  examined.  This 
hypothesis  works  better  in  the  trough  where  turbulent  bottom  boundary  layer  processes  are 
more  dominant  than  over  the  bar,  where  breaking-wave  induced  turbulence  generated  at  the 
surface  modifies  the  profile. 

A  relationship  between  near  bottom  vertical  velocity  profiles  and  bottom  roughness 
was  found.  While  the  relatively  small  amount  of  data  used  (three  days)  are  not  enough  to 
establish  a  definitive  conclusion  from  the  available  data,  it  can  be  concluded  that  the  bed 
shear  stress  coefficient  is  directly  proportional  relationship  to  the  bottom  roughness  and 
therefore  is  an  important  parameter  to  characterize  the  bottom  boundary  layer.  Surprisingly, 
only  poor  correlation  of  and  was  found  with  wave  parameters. 

The  influence  of  the  wind  is  parameterized  by  an  exponential  approximation  that 
works  well  in  50%  of  the  analyzed  profiles.  This  parameterization  assumes  that  the  residual 
between  the  velocity  corrected  using  the  surface  elevation  pdf  and  the  measured  data  is  only 
due  to  the  wind  and  does  not  consider  the  alongshore  component  of  mass  transport  velocity 
that  is  assumed  small. 

These  results  should  be  treated  carefully  since  there  were  obtained  from  a  small 
number  of  days  over  a  barred  beach  with  strong  longshore  currents  and  cannot  be 
generalized  without  further  studies. 
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APPENDIX  A.  FIGURES 
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Figure  2.  Bathymetry  and  equipment  positioning  during  DUCK94  experiment  (11  Oct.) 
Sled  transect  line  was  940  m  alongshore 


Profile  Line  Numbers 


Depth  (m) 


Angle  0  T  (sec 


Figure  4  Climatology  for  the  three  days  considered  (10,  11  and  12  Oct  ).  Current 
measured  in  the  middle  of  the  trough.  H  is  the  significant  wave  height  and  T  is  the 
of  peak  frequency 


DUCK  94  -  Longshore  current  -  Oct  11  time  15:44-  xposition  =  157.3 


Figure  6.  Example  of  vertical  profile  (the  upper  two  current  meters  sometimes  were  out 
of  the  water). 
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DUCK94  -  Bottom  roughness  -  Oct12,  Line  10 


Figure  8.  Variation  of  the  Rtiis  Bottom  Roughness  with  the  Cross-Shore  Distance 


DUCK94  -  Cf  vs  Um/v 
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DLICK94  -  Cf  vs  Bottom  roughness 


Figure  1 1  Bed  Shear  Stress  versus  Bottom  Roughness.  The  Line  Represents  a  Linear 
Regression  Based  on  All  Points  but  the  Outlier. 


z  ^ 

Figure  13  —versus  —  (Outlier  not  Included). 
h  A 


Eddy  viscosity  (cm2/s) 


DUCK94  -  Eddy  Vis.  vs  X-distance 


X-distance  (m) 


Figure  15  versus  Cross-Shore  Distance 
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DUCK94  -  Cf  vs  Gamma 


Figure  16.  Relationship  Between and  y 


Eddy  Viscosity  (cm2/s) 


DUCK94  -  Eddy  Viscosity  vs  Gamma 


Gamma 


Figure  17  Relationship  Between  and  y 
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DUCK  94  -  Longshore  current  Oct  12  time  11:2  xposition  =  251.5 


Figure  18.  Example  of  wind  contribution  to  the  longshore  current:  exponential 
approximation  to  the  residual. 
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X-Pos. 

(m) 


(cm/s) 


(cm  ^/s) 


101 

7:33 

291 

3.80 

1.11 

0.29 

0.03 

2.70 

0.002 

68.3 

0.98 

102 

9:06 

261 

3.24 

1.12 

0.35 

0.01 

1.85 

0.001 

39.9 

0.96 

103 

10:10 

231 

3.26 

1.15 

0.35 

0.50 

5.62 

0.005 

122.3 

0.97 

104 

12:22 

201 

3.45 

0.89 

0.26 

0.67 

7.18 

0.005 

165.1 

0.99 

105 

14:06 

185 

3.18 

0.75 

0.24 

2.57 

7.09 

0.011 

150.4 

0.99 

106 

15:36 

170 

2.57 

0.63 

0.25 

1.37 

3.18 

0.009 

54.5 

0.99 

107 

16:57 

156 

1.98 

0.55 

0.28 

0.82 

2.13 

0.014 

28.1 

0.99 

111 

7:53 

292 

3.66 

1.30 

0.36 

0.28 

5.92 

0.004 

144.3 

0.96 

112 

9:29 

270 

3.16 

1.27 

0.40 

0.23 

5.02 

105.9 

0.96 

113 

10:45 

244 

3.16 

1.12 

0.35 

1.10 

5.85 

0.006 

127.4 

0.97 

114 

12:01 

230 

3.75 

0.98 

0.26 

19.98 

14.28 

0.047 

357.4 

0.91 

115 

13:17 

187 

3.24 

0.89 

0.27 

1.23 

5.29 

0.007 

114.4 

0.99 

116 

14:30 

171 

3.05 

0.84 

0.28 

3.05 

4.05 

0.011 

82.3 

0.99 

117 

15:44 

•  157 

2.50 

0.76 

0.30 

0.93 

3.01 

0.009 

50.1 

0.99 

118 

17:03 

146 

1.63 

0.66 

0.40 

3.23 

2.25 

0.048 

24.4 

0.89 

121 

7:44 

298 

3.51 

1.23 

0.35 

0.03 

1.70 

0.002 

39.9 

0.97 

122 

9:27 

273 

2.93 

1.30 

0.44 

0.03 

2.44 

0.003 

47.7 

0.93 

123 

11:02  ^ 

252 

2.98 

1.21 

0.41 

0.43 

4.23 

0.005 

84.1 

0.99 

124 

12:26 

225 

3.54 

1.04 

0.29 

6.93 

8.11 

0.018 

191.5 

0.99 

125 

13:26 

210 

3.51 

0.97 

0.28 

2.37 

6.17 

0.010 

144.5 

0.99 

126 

15:13 

188 

2.99 

0.90 

0.30 

1.83 

4.72 

0.009 

94.0 

0.99 

127 

16:22 

172 

2.74 

0.84 

0.31 

7.05 

3.26 

0.020 

59.6 

0.98 
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Run 

Time 

Cross-Shore 
Position  (m) 

Wind  speed 
(m/§ 

Wind  Cuixent 
angle(  ) 

101 

7:33 

291 

12.2 

23 

102 

9:06 

261 

12.1 

16 

103 

10:10 

231 

12.9 

19 

104 

12:22 

201 

11.7 

15 

105 

14:06 

185 

11.4 

11 

106 

15:36 

170 

11.6 

9 

107 

16:57 

156 

11.2 

12 

111 

7:53 

292 

12.3 

44 

112 

9:29 

270 

11.8 

40 

113 

10:45 

244 

10.9 

39 

114 

12:01 

230 

11.0 

39 

115 

13:17 

187 

11.1 

35 

116 

14:30 

171 

11.3 

40 

117 

15:44 

157 

10.6 

41 

118 

17:03 

146 

11.0 

41 

121 

7:44 

298 

14.1 

56 

122 

9:27 

273 

14.3 

56 

123 

11:02 

252 

13.5 

55 

124 

'  12:26 

225 

14.1 

57 

125 

13:26 

210 

13.9 

59 

126 

15:13  : 

188 

13.9 

59 

127 

16:22 

172 

13.8 

58 

Table  2.  Wind  Data 
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